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^ : Abstract 
in 

We show that the one-loop QCD correction to deeply-vhtual Compton scat- 

^ . tering can be factorized into the finite perturbative contributions and the 

jy-^ , collinearly-divergent terms, which correspond to the matrix elements of the 

I off'- forward parton distributions. As a by-product, we obtain the next-to- 

Q I leading order coefficient functions in the generalized operator product expan- 

■ sion of two vector currents. 
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In searching for ways to measure the amount of the nucleon spin carried by quark orbital 
angular momentum, one of us introduced deeply-virtual compton scattering (DVCS) as a 
probe to a novel class of "off-forward" parton distributions (OFPDs) [|I|. DVCS is a process 
in which a highly virtual photon (with virtual mass y> ^qcd) scatters on a nucleon 
target (polarized or unpolarized), producing an exclusive final state consisting of a high- 
energy real photon plus a slightly recoiled nucleon. With the virtual photon in the Bjorken 
limit, a quantum chromodynamic (QCD) analysis shows that the scattering is dominated 
by the simple mechanism in which a quark (antiquark) in the initial nucleon absorbs the 
virtual photon, immediately radiates a real one, and falls back to form the recoiled nucleon. 

Several interesting theoretical papers have since appeared in the literature, which made 
further studies of the DVCS process f^-^] and the OFPDs A process closely-related 

to DVCS, Compton scattering with two- virtual photons, have been studied before from 
theoretical interests The OFPDs have also been recognized in a number of theo- 



retical studies in the past |T^,[TB|. Moreover, the QCD evolution equations for the OFPDs 
PJ3|J^,P|,PH1 are closely related to the evolution of meson wave functions and light-ray or 



string operators However, physical significance and actual measurement of OFPDs 

have not been thoroughly explored in the literature. It was pointed out in Ref. [0] that 
the OFPDs can appear in general types of hard diffractive processes. Recent works in Ref. 
||T5Hl8[ have indeed found their use in diffractive meson production |19,20] and Z° or muon 
pair production pT],E2|. 



In this paper, we study the one-loop correction to DVCS in QCD. Here, for simplicity, we 
consider only the parton helicity-independent amplitude. [The helicity- dependent case will 
be published separately, together with a more detailed analysis of the present calculation 
23| .] Our result shows that the infrared divergences in the one-loop diagrams can be entirely 



factorized into the nonperturbative OFPDs. This result is not immediately obvious from an 
analogy with a general virtual Compton scattering, because in DVCS the final state photon 
is real and new infrared divergences can potentially arise from the appearance of a new light- 
like momentum. It turns out, however, that a factorization theorem does hold at one- loop 
level with the final state photon real and structureless. We expect the factorization theorem 
to be true to all orders in perturbation theory, allowing DVCS to be studied in perturbative 
QCD just like other gold-plated examples such as deep-inelastic scattering. The finite part 
of our one-loop results, together with the two-loop evolution of OFPDs in the MS scheme, 
provides the necessary ingredients for calculating DVCS at the next-to-leading order. 

To carry out our study in a more general setting, we actually consider the non-forward, 
unequal-mass virtual Compton scattering |5|, p!0| , pJ]| . We call the incoming (out-going) virtual 



photon momentum [q'^ = q^ — A'^), and the incoming (out-going) nucleon momentum 
pM (^p'M = _|_ A^). We introduce the "average" photon momentum q^ = {q + q')^ /2 and 
"average" proton momentum = {P + P'Y /2. The Compton amplitude is defined as. 



T^"" = i d'^ze-'^-' ( P' 



P) ■ (1) 



We want to study the simplification of the amplitude in the Bjorken limit: = —(f oo, 
P ■ Q —>■ oo, and the ratio of the two staying finite. 

For convenience, we choose two light-like momenta and so that = n'^ = and 
p ■ n = 1. We let the average photon and nucleon momenta to be parallel to and n^, so 
that. 
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where = — and M is the mass of the nucleoli. The Bjorken xb is Q"^/ (2P • q), 
where corrections of order 1/Q^ are ignored. Clearly, any four-momentum can be expanded 



in terms of 



and two unit vectors in the transverse directions. Large scalars in any 



Feynman diagram come from the product of the p component of a four-vector and the n 
component of q. Therefore, in the leading order, we can safely ignore other components of 
an external momentum. For instance, the momenta of the initial and final state nucleons 
can be approximated as (1 + and (1 — C,)p^, respectively, where ^ is defined form the 
expansion, 



A^' = -2i p^" + 



(3) 



and is constrained to [0,1] by our choice of coordinates. [Note that the definition of ^ differs 
from that in Ref. by a factor of 2.] Moreover, the momenta of initial and final partons 
participating in a hard subprocess can effectively be taken as (x -|- ^)p^ and (x — ClP^i 
respectively, where — 1 < a; < 1. The longitudinal momentum fractions of the initial and 
final nucleons carried by partons are (x + ^)/(l + ^) and (x — ^)/(l — ^), respectively. 

We consider only the symmetric part of T^'^, which is insensitive to helicities of partons 
in a nucleon target. Like the forward Compton scattering, there are two leading-order 
Lorentz- invariant amplitudes [H, 



q'lJ-p . q' 
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where Tj are functions of xb, ^, t = A^, and Q^. Physically, Ti represents the amplitude for 
the transversely-polarized photon scattering, and the amplitude for the longitudinally- 
polarized scattering. Since does not contribute to DVCS, we ignore it in the remainder 
of the paper. 

A factorization theorem is believed to exist for the general virtual Compton amplitudes 
defined above. For instance. 



/I dx — / 
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where a labels different parton species. Fa is the helicity independent, leading-twist off- 
forward parton distributions as defined in Ref. For quarks, one has. 



Fa=,ix) 



dX 
— ( 
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and similarly for gluons, Fg{x). Cia is the coefficient function calculable in perturbative 
QCD. At the leading order in one has 
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where Ca is the electric charge. The second term in the bracket represents a crossing con- 
tribution which is always present in Compton processes. To simplify the presentation, we 
will not include it explicitly in the following formulas. For ^ = t = 0, we go back to the 
well-known forward virtual Compton scattering. 

The factorization formula works for both xb > ^ and < 1 regions. In fact, it is 
defined throughout the complex plane of a;^- The amplitude is analytic at \xb\ > 1, where 
there is no on-shell propagation. Therefore, one can expand Cia in Taylor series at xb = oo, 

cU-,h= E CTMQ'))-^. (8) 

XB XB n=2,4... 1=0,2... B 

The x-integration in Eq. can now be done using 

dxx''-^Faix,i,t) = aJi,t) = -n^'^n^^---ni'-iP'\ipiD ■■■iD 'y''"MP). (9) 
-1 2 

Therefore the i = terms in the amplitude Ti are just the result of the usual twist-2 
operators in the operator product expansion (OPE) for forward virtual Compton scattering. 
As for the i terms, one can translate the ^ factors into total derivatives on the twist-2 
operators, 

n^,...n^„+AP'\^d^'"+\..id^-+'ijz D ---iD j^"ij\P) = 2{20''an{^,t) . (10) 

Therefore the factorization formula in Eq. (^ reflects a general OPE in terms of twist-2 
operators with arbitrary numbers of total derivatives. C^^lasiQ"^)) is the coefficient function 



relevant for deep-inelastic scattering |2J]. A general virtual Compton process also requires 
the i 7^ terms. 

Deeply virtual Compton scattering is a physical process with the kinematic requirement 
that the final photon be on its mass shell. This corresponds to the region x^ = ^ < 1, 
where the OPE form of the amplitude is not particularly illuminating. In particular, the 
Compton amplitudes now have both real and imaginary parts. The analytic continuation 
to the \xb\ < 1 region is made by approaching the real axis from the lower half plane, and 
hence xb has a small negative imaginary part. The central question we would like to address 
is this: Does the factorization theorem still hold when = ^? One can, of course, study 



this question using the general factorization techniques [|25|. While we believe the answer is 
affirmative, we present a one-loop calculation to support this. 

We have calculated the general one-loop virtual Compton scattering with on-shell quarks 
and gluons. We present the quark result first. The initial and final momenta of the quark 
target are taken to be = (x + ^)p^ and P'^ = (x — ^)p^, respectively. For convenience, 
we replace the spinor-space matrix element u{P')Tu{P) by a trace Tr[|(r]/2. Although they 
are not identical, this replacement only affects the interpretation of the result. We call the 
resulting Compton amplitude, t^^ , with Lorentz-scalar amplitudes ti. At the leading order, 
one has, 

ha = . (11) 

X — Xb + 26 

A practical aspect of parton calculation lies in the fact that the coefficient function Cia(x, ^) 
can be obtained from the infrared subtracted parton scattering amplitude. 
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We use dimensional regularization to regularize the infrared divergences present in the 
loop calculations = 4 + e, e > 0). We use the Lehmann-Symanzik-Zimmermann reduction 
formula to extract on-shell matrix elements from the corresponding Green's functions. The 
calculation is done in the region where > 1 so that there is no imaginary part. Sum- 
ming up the standard vertex, self-energy, and box diagrams, we find the following divergent 
contribution, 



.pole 
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An examination shows that the above expression is proportional to the perturbative ma- 
trix elements of the OFPDs calculated in dimensional regularization at one- loop order 0. 
Therefore, the effect of this term is already taken into account in the general factorization 
formula with a leading-order coefficient function C^^ and the nonperturbative OFPDs. In 
the limit xb —>■ C,, the factorization is not affected. 
We also find the finite part of the one- loop result. 
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(x^ - - LO^x^) (x^-e^) \2UJ^ l-UJ^X^ 

'{l-ujxf + 2{ujx-uj'^i'^) 
2uj{x'^ - ~ ujx) 



3 2x-u{x^ + e 
iog(l — uox) — -- 



2{x'^ -i'^){l-ujx) 



log(l - uji) 
log(l-a;x)| , (13) 



where we have introduced uj = 1/xb- The above expression is manifestly finite in the DVCS 
limit, uC, ^ 1, so the one-loop factorization theorem for DVCS holds for quark scattering. 
When ^ = 0, we recover the coefficient function for the forward Compton scattering pertinent 
to deep-inelastic scattering |^ . The coefficient function coupling with the off- forward quark 
distributions in the general virtual Compton scattering at order is. 



c|,(x,o = xtff^(x,e,c^ = i 



(14) 



Next, we consider virtual Compton scattering on an on-shell gluon. Again for conve- 
nience, we replace the polarization product e'^*e'^ by {—g^^ + p^n^ + p'^n^^)/{2 + e). We 
neglect the last two terms because of the color gauge invariance. At one-loop level, there 
are six quark-box diagrams, three of which have identical results as the other three. Two 
of the three inequivalent diagrams are related to each other by crossing symmetry; and the 
remaining one is itself crossing-symmetric. The infrared divergent part of the amplitude is. 
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where Tj? = 1/2. It is not immediately clear that all the terms above can be absorbed into 
the renormalization mixing of the quark distributions with the gluon. However, notice that 
the helicity- independent gluon distribution Fg{x) is antisymmetric in x — >■ —x. Therefore, 
when convoluted with Fg{x) only the x-symmetric part of the amplitude contributes. And 
the x-symmetric tig^'^ is identical to the mixing coefficient in a renormalized quark density 

i- 

The remaining finite part is considered as the gluon contribution to the Compton scat- 
tering in the MS scheme. Our calculation yields. 



t 



finite 




1 
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(17) 



The above expression is finite in the limit of cj^ ^ 1, so the factorization for DVCS on 
a gluon target holds at one-loop order. For ^ = 0, our result agrees with the coefficient 
function calculated in ||2^. Using the definition of the off-forward gluon distribution Fg{x), 



we found the coefficient function in the general factorization formula at the as order. 



cUx,o 



X 



x^ 



finite / 



{x,i,uj = 1) . 



(18) 



In a recent paper, Miiller |26| studied the constraints of the spacetime conformal sym- 
metry on the form of the OPE for two vector currents. He found that at the next-to-leading 
order, the OPE can be expressed in terms of the coefficient functions in the forward scat- 
tering (^ = 0) and the improved conformal-covariant operators at the leading order. While 
the result can be checked for the non-singlet case using the two-loop anomalous dimensions 
in the literature |27], such anomalous dimensions for the singlet case are not yet available. 



Clearly, they are needed if one is interested in carrying the analysis for DVCS to the full 
next-to-leading order. 
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